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Abstract 

We present a chiral model for the J = 0, / = 1/2, elastic Ktt amplitude, suited to be employed in 
K^TT^TT^ data analyses and valid between threshold and 1.5 GeV. Although not as precise as 
other versions available in the literature, it is rather simple and incorporates the essential physics 
in this energy domain. In the case of the iC-matrix approximation, the model allows the pole 
structure of the Kir amplitude to be understood by solving a quadratic equation in s. We show 
that the solutions to this equation can be well approximated by polynomials of masses and coupling 
constants. This analytic structure allows a clear understanding why, depending on the values of 
one of the coupling constants, one may have one or two physical poles. The model yields a pole, 
associated with the k, at y/s = (0.75 — iO.24) GeV. 
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I. INTRODUCTION 



Since the E791 experiment [1], heavy-meson decays have systematically produced solid 
evidence in favour of low-energy scalar resonances [2-4]. The quahty of these results renewed 
interest in the problem and motivated an effort aimed at determining the positions of their 
poles in the complex energy plane. In the case of tttt resonances, the pole position of the 
a can be determined from scattering data by means of the Roy equation, and the rather 
precise value = (0.44ll8^ - i 0.272li2.5) MeV is available[5]. 

The situation of scalar resonances in the K'k system is much less certain, since our knowl- 
edge of the 5"- wave 1 — 1/2 amplitude is based on just two experiments [6, 7], which include 
phase shifts only up to ~ 80°. These data sets have been carefully analyzed in recent years, 
and determinations of both pole positions[8] and scattering lengths[9] became available. 
Nevertheless, the need of more empirical knowledge is still urgent and, in principle, informa- 
tion from the decay K~t:^t:^ could be useful in either constraining or complementing 
Ktt scattering data. 

Analyses of decays normally rely on trial functions written in terms of Breit-Wigner 
expressions, which are at odds with chiral symmetry. As the symmetry is very important at 
low-energies, this kind of procedure has been criticized and alternatives were proposed[10, 
11]. In the work by 011er[ll], data were refitted with the help of a chiral amplitude, with 
a remarkable decrease in x^. In a recent paper [12], we have discussed the main features 
of the decay K~t:^t:^ at low energies, including explicit descriptions of both the 

primary weak vertex and final state interactions, based on a unitarized Kn amplitude. This 
elastic amplitude was obtained by the iteration of a kernel by means of a simplified Bethe- 
Salpeter equation[13]. The (J,/) = (0,1/2) component of the kernel, which interests us 
here, was derived from effective lagrangians and contains a leading 0{q^) contact term[14], 
supplemented by an explicit resonance exchange[15], corresponding to 0{q^) corrections. 
This model is illustrated in fig.l. 

Our ultimate goal is to produce a tool to be directly employed in analyses of raw data 
and therefore we try to avoid, as much as possible, the use of long and involved expressions. 
We thus neglect contributions from resonance exchanges m.t— and u— channels, since they 
just give rise to small backgrounds [12]. In the present work we follow this approach and 
show that simplifications are possible in the chiral kernel which preserve its essential physical 
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FIG. 1: (Color online) Tree-level txK amplitude; dashed and full lines represent respectively pseu- 
doscalar mesons and scalar resonances. 
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II. Ktx amplitude 



Our unitarized (J, /) = (0, 1/2) amplitude for the process n K n K has been discussed 
in detail in ref.[12] and here we just summarize its main features. It is written as 



l\s)/D{s) 



D{s) = [ml-s + j\s) R^Js)] 



lOTT 



(1) 



where: 



- s is the usual Mandelstam variable and p{s) = a/1-2 {Ml + M^)/s+ {M^-M^y/s'^ ; 

- mji is the parameter present in the chiral lagrangian, called nominal resonance mass; 

- R^^^{s) is the function describing off-shell effects in the two-meson propagator, given by 

R^^^is) = -^[Lis)~Liml)]/lQn\ 
^L{s) = p(s)log[(l-(7)/(l + a)]-2 

+ [(Mi-M2)/s] log(M^/MO], 



a = ^/\s-{MK + M^)^\/\s-iMK-M^y\ ; 

-^1/2 ("^i?) = t)y construction and therefore the phase shift is 7r/2 at s = m|j ; 
7^(s) is the function which incorporates chiral dynamics, given by 

7^(.) = {(l/F^) [(l-3/(s)/8). 

+ {(3/F^) [q {s-M^M],) 
+ Cm (4M|, + 5M3 /6]2} 



(2) 



R 



(3) 



- the labels L and R in the curly brackets denote contributions from the leading 0{q^) 
contact term and the 0{q^) resonance correction; 

- F is the pseudoscalar decay constant and the parameters and are the resonance 
couplings defined in ref. [15]. 

This model for the Kt^ amplitude depends on six parameters, three of which are well 



known and given by = 0.1396 GeV, Mk = 0.4937 Gev and F = y/F\F^ = 0.103 GeV. 
The other three, namely m^, q and need to be taken from experiment. In ref. [15] one 
finds \cm\) = (0.032, 0.042) GeV, obtained from the process rjn. We adopt these 
values provisionally and, at the end, suggest our own choice. 




s(GeV^) 



FIG. 2: (Color online) Modulus of the elastic amplitude: the full curve represents eq.(3), the 
dashed curve L is obtained by making c^^ — Cm — 0, and those labeled and are resonance 
contributions proportional to these parameters. 



The infiuence of the various dynamical mechanisms over \T_^^J is displayed in Fig. 2. The 
full curve is based on eq.(3), whereas the dashed ones describe partial contributions: L 
corresponds to q = Cm = and those labeled Cd and Cm arise from just the terms proportional 
to these parameters. The curve Cm has a typical Breit-Wigner shape, but that labeled Cd 
is much wider. Nevertheless, both of them are constrained to be small at threshold, since 



4 



the resonance corresponds to a 0{q'^) correction. The leading 0(0^) contact term clearly 
dominates at low-energies and, as expected, full and resonance curves coincide at s = m|j, 
since the resonance nominal mass is independent of coupling parameters. By construction, 
T is purely imaginary at this point. 



III. POLES 



The poles of T^^^j determined by the condition D{s) — 0, can be easily found out by 
numerical methods. Before doing this, however, it is instructive to discuss the mathematical 
structure of the problem, by means of a simplified version of the amplitude, in which: 

- the function R^^^ is neglected, which corresponds to the X— matrix approximation; 

- the pion mass is neglected, which corresponds to the SU{2) limit. 
The poles can then be found by solving the complex quartic equation 



5/8- 



F2 



+ [-{5ml + 7M|.)/8 



Cd 



Mi 



+[{ml + 3M|,)/8 + 3(crf - 2cJ^fMl/ F'']M^^ s 
-^m\M%/S = . 



(4) 



Around physical poles, the quantity M|^/|sp is small and one obtains the quadratic 
equation 



As^ + Bs + C^O, 
A = [5/8 - Scl/F'] , 
B = [-{5ml + 7Ml)/l 
C ^[7M]^/8-il67rF^]ml 



(5) 



Q(9Q-4cJ^ + il67rF2] 



The parameter A plays an important role in this problem and receives contributions from 
both the leading contact term and the resonance. In particular, the condition ^4 = yields 
Cd — F-\Jb/2A — 0.047 GeV, which is jut 50% larger than accepted empirical values. In the 
sequence, we discuss the behavior of the solutions of this equation as functions of A, in the 
interval 5/8 < A < 0. Exact solutions are available in its two extremes: 



- A — 5/8 ^ Cd — : in this case, the resonance R corresponds to a bound state in the real 
axis, which does not couple with the Ktt system, and one has 



s+(5/8) = m% and S-(5/8) = [7M^/5 - i USnF^/b] ; 

^4 = : the dynamical equation is no longer quadratic and its single solution reads 

[7M|:/5-i 1287rF75] 



._(0) 



1 + 



7M' 



Ml 



^-8Q(9Q-4c^)9^-il287r^ /m% 



(6) 



(7) 



We note that these two solutions already show a prominent feature of the problem, 
namely the stability of the solution s_ (^4) in the whole interval considered. As is a large 
parameter, s_(0) = S-(5/8) + 0(l/m|). 

In the general case, the solutions of eq.(5) will have the form 



s± ^[-B±^/D ]/2A ^ D ^ - AAC . 



(8) 



The square root prevents algebraic simplification of results. However, at the point A = 
—'SlB/2m\ ~ 5/16, one has Q^D — and an approximate solution can be obtained for 
. By imposing a quadratic polynomial in A to interpolate this function at A = 5/8, 
—'RB/2m\ and , one finds 
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(9) 



(10) 



In tables I and II we display figures derived from eqs.(4), (5) and (9-10), for a wide sample 
of values for Cd/F. We recall that empirical estimates of this quantity lie in a narrow band 
around Cd/F ~ 0.3 and the condition A — corresponds to q/F = 0.456. As far as accuracy 



6 



is concerned, one learns that predictions from the quartic and quadratic equations deviate 
very little. The approximate algebraic solutions are less precise but, nevertheless, describe 
well the qualitative features of the results in the whole range considered and are reasonably 
precise in the region of empirical interest. The accuracy of cqs.(9-10) could be improved 
by keeping more terms in series expansions, but this would just yield larger expressions, 
without further conceptual gains. 

TABLE I: Poles in GeV, of the amplitude T^^^, for niR = 1.2 GeV and = 0.042 GeV. 





quartic 


quadratic 


analytic 


simplified 


Cd/F 


eq.(4) 


eq.(5) 


eq.(9) 


eq.(ll) 





1.2000 


1.2000 


1.2000 


1.2000 


0.1 


1.2351 -z 0.0269 


1.2357 -z 0.0272 


1.2394 -z 0.0269 


1.2300 -z 0.0175 


0.2 


1.3135 -i 0.0898 


1.3155 -i 0.0918 


1.3257 -z 0.0874 


1.3371 -z 0.0758 


0.3 


1.5346 -i 0.2501 


1.5346 -i 0.2500 


1.5513 -z 0.2245 


1.6050 -z 0.2022 


0.4 


2.4694 -i 0.6287 


2.4644 -i 0.6327 


2.4566 -z 0.6408 


2.5521 -z 0.5534 


0.5 


0.8660 + i 2.7361 


0.8586 + i 2.7414 


1.0882 + z 2.8707 


0.9040 + z 2.8315 



TABLE II: Poles ^sT, in GeV, of the amplitude T^^^, for tur = 1.2 GeV and Cm = 0.042 GeV. 



Cd/F 


quartic 

(x-i.(-]) 


quadratic 


analytic 

c'q.(lO) 





0.7908 -z 0.5294 


0.7908 -z 0.5294 


0.7908 -z 0.5294 


0.1 


0.8064 -z 0.5109 


0.8014 -z 0.5173 


0.8001 -z 0.5242 


0.2 


0.8458 -z 0.4819 


0.8390 -z 0.4879 


0.8206 -z 0.4848 


0.3 


0.8904 -z 0.4085 


0.8908 -z 0.4126 


0.8586 -z 0.4215 


0.4 


0.8756 -z 0.3228 


0.8837 -z 0.3165 


0.9189 -z 0.3391 


0.5 


0.8440 - z 0.2736 


0.8576 -z 0.2575 


1.0041 -z 0.2459 



Conceptually, the behaviours of y^s+ and ^/sZ shown in the tables are strikingly different. 
The latter is a rather slow- varying function, whereas the former changes rapidly and even 
has the sign of the imaginary part reversed in the last row. In order to understand this 
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behavior of ^/sljl, we keep just the leading term in eq.(9) and find 

In spite of its simphcity, this result yields quite reasonable predictions, as indicated in the last 
column of table I. Moreover, it shows clearly that the major features of s+ are determined 
by the factor A in the denominator. In particular, it explains the change in the sign of 
the imaginary part of ^/slfT observed in the bottom line of the table, since A vanishes at 
Cd/F = 0.456. As this factor does not occur in s_, table II is much more monotonic. 
The following scenario is supported by eqs.(9-10): 

- in case the resonance R is absent, one has just the pole y^sT, which is due to the leading 
contact interaction; 

- if the resonance is present, but its couplings to mesons are not turned on {cd — Cm — 0), 
one has the pole y^sT and a bound state in the real axis at s = m|j[26]; 

- when the parameters and c„i are turned on and the resonance R couples to ^/sZ, both the 
mass and width of y/s+ increase monotonically, driven by the factor A in the denominator; 

- as a consequence, for realistic values of q, one has necessarily ^^/s~^ > tur and the pole 
on the complex plane has to stand on the right of the point at which the phase shift passes 
through 7r/2; 

- the pole Y^s+ blows up at the critical value c^/F = •\/5/24 and, beyond this point, just 

^JsZ is present; 

Although the couplings of the resonance R to mesons is implemented by both q and 
Cmi the former is much more important than the latter, for it is incorporated into the 
parameter A. The term proportional to Cm is less relevant, because it contains meson 
masses and vanishes in the chiral limit. For instance, if one chooses Cd/F — 0.3 and sets 
Cm — 0, the solutions of the quadratic equation become = 1.3631 — i 0.1973 GeV and 
y/sT — 0.9971 — i 0.4835 GeV. This means that gentle variations of around empirical 
averages would have little influence over pole positions. 

In figs. 3 and 4, predictions from the full model for E = ^/s, given by eq.(l), are com- 
pared with results from K-matrix and the quadratic approximations. Quartic and quadratic 
approximations cannot be distinguished visually and the former is not shown. Inspecting 
these figures, one learns that the inclusion of the pion mass is not numerically important, 
but off-shell effects in the two-meson propagator do infiuence the positions of the poles and 
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tend to decrease both masses and widths. Nevertheless, it does not alter the qualitative 
features of the scenario discussed above. 
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FIG. 3: (Color online) Real (full lines) and imaginary (dashed lines) components of the function 
= ^JsZ; the full curve comes for eq(l), whereas those labelled K and q correspond to the 
K-matrix and quadratic approximation. 




Our results suggest that the pole ^Js^ can be identified with the K^{1A?)Q). The fitting of 
both the mass and width of this state supplies two constraints for the resonance parameters. 
As Cm has little influence over numerical results, we flx it at the value Cm = 0.042 GeV and 
choose tur and q so that ^ = [(1.414 ± 0.006) - i (0.145 ± 0.010)] GeV [27]. We then 
flnd Cd/F = 0.2705 ± 0.0078 , niR = 1.1865 ± 0.079 GeV and ^ = (0.7505 ± 0.0010) - 
i (0.2363 ± 0.0023) GeV. It is interesting to compare our results with those from the much 
more complete analysis by Jamin, Oiler and Pich[8]. Although there are no data in the 
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FIG. 4: (Color online) Real (full lines) and imaginary (dashed lines) components of the function 
= y^s^; the full curve comes for eq(l), whereas those labelled K and q correspond to the 
K-matrix and quadratic approximation. 

(J,/) = (0,1/2) channel for phase shifts around 7r/2, parametrizations produced in their 
paper suggest that rriR ^ 1.32 GeV and our value misses this point by 10%. On the other 
hand, by fixing the mass of the Kq and imposing mji = 1.32 GeV, we find q/F = 0.1640, 
= 0.7419 — i 0.2541 GeV, and Tk* = 0.1035 GeV, which is too small. In both cases, 
the values of Cd/F are below those of ref. [15]. The difficulty of fitting simultaneously rriR 
and Tk* may be associated with the fact that we did not include inelasticities present in the 
region s > 1 GeV. 

For the sake of completeness, in table III, we compare our result, ^/sZ = 0.751 — 
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i 0.236 GeV, with the positions of the k pole obtained in previous works. 



TABLE III: Pole position of the k, in GeV. 



year 


ref. 


K 


1986 


[16] 


0.727 0.263 


1997 


[17] 


0.905|g 0.222^^5 


1998 


[18] 


0.911 -i 0.158 


1999 


[19] 


0.779 -i 0.330 


2000 


[8] 


0.708 0.305 


2002 


[2] 


0.721 ± 19 ± 1-1 - / 0.292 ± 21 ± 44 


2003 


[20] 


0.750^ 0.342 ±60 


2004 


[21] 


0.750 ± 18 -i 0.226 ±11 


2006 


[22] 


0.694 ± 53 -i 0.303 ±30 


2006 


[23] 


0.658± 13-z0.279± 12 


2006 


[24] 


0.841 ± 23^1 0.306 ±26^J 


2008 


[25] 


0.772 -i 0.281 



IV. SUMMARY 

In this work we have discussed a model for the Kn amphtude in the (J,/) = (0, 1/2) 
channel, suited to energies up to 1.5 GeV. It is aimed at being used in data analyses of 
processes such as D'^ — > K'Tr'^Tr'^ and given by 

„2 
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1/2 



(m|-s) -i7V/(167r) 



(12) 



+ s 



Cd 



F2, 



^ (, _ Mi) + i|^M- 



{m\ - s) 



It represents a compromise between simplicity and the essential phenomena of this channel. 
The factor (s — M|-) in 7^ is the Adler zero, which arises naturally in the framework of chiral 
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symmetry. It allows automatically for one complex pole, identified with the n. Depending 
on the values obtained for the free parameters niR, ca and c^, another pole can be present, 
associated with the Kq. 
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